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We study the single transverse spin asymmetry of prompt photon production in high energy 
proton-proton scattering. We include the contributions from both the direct and fragmentation 
photons. While the asymmetry for direct photon production receives only the Sivers type of con- 
tribution, the asymmetry for fragmentation photons receives both the Sivers and Collins types of 
contributions. We make a model calculation for quark-to-photon Collins function, which is then used 
to estimate the Collins asymmetry for fragmentation photons. We find that the Collins asymmetry 
^■s^- for fragmentation photons is very small, thus the single transverse spin asymmetry of prompt pho- 

ton production is mainly coming from the Sivers asymmetry in direct and fragmentation photons. 

1 We make predictions for the prompt photon spin asymmetry at RHIC energy, and emphasize the 

. importance of such a measurement. 

PACS numbers: 24.85.+P, 12.38.Bx, 12.39.St, 13.88.+e 
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I. INTRODUCTION 

-P 

Single spin asymmetries (SSAs) in transversely polarized proton-proton collisions have provided essential informa- 
tion on the internal partonic structure of the proton, particularly the parton's transverse motion in the transversely 
polarized proton Two different yet related QCD factorization formalisms have been proposed to describe the 
,-P ■ observed asymmetries: the transverse momentum dependent (TMD) factorization and the collinear twist-three 
l— ~" \ factorization approaches (tHUJ- 

■ For processes such as semi-inclusive hadron production in lepton-proton deep inelastic scattering (SIDIS) lp^ — > 
" I'hX which are characterized by both the photon virtuality Q 2 and hadron transverse momentum Phi. such that 
Q ^> Ph± ~ Aqcd, one describes the SSAs in the TMD factorization formalism. In this approach the transverse 
spin effects are associated with naive time-reversal-odd TMDs which represent helicity flip cut quark target scattering 
amplitudes with a non-trivial color phase flij . Two well-known TMDs are the quark Sivers function (l5^ and Collins 
function [l6j], which describe the so-called sin(</>/j — S ) and sin^ + (f> s ) modulations in SIDIS on transversely polarized 
OO ] target, respectively. Because of the different angular modulations in the cross section, one can separate Sivers from 
Collins effect in SIDIS and thus extract them independently from the experimental data [l7l - [l9| . On the other 
hand, for single inclusive hadron production in proton-proton scattering p^p — > hX where there is a single hard 
scale given by the hadron's transverse momentum, Phi. 3> Aqcd, one can describe the SSAs in the collinear twist- 
three factorization approach in terms of either the twist-three quark-gluon correlation functions in the transversely 
polarized proton [20l l2l| , or the twist-three fragmentation functions in the hadronization process [lH ] . In the twist- 
three formalism we refer to the former contribution as Sivers effect, and the latter one as Collins effect, since they 
represent the collinear version of these two effects (based on the operator definitions relating the first fey-moments 
of the Sivers and Collins functions in the collinear twist-three approach Q). While the most abundant experimental 
data exist (22l . l23j on transverse spin effects in SSAs of single inclusive hadron production in proton-proton collisions, 
disentangling Sivers and Collins contributions presents a significant experimental challenge, thus the true origin (the 
relative contributions from these two effects) for the inclusive hadron production still remains elusive [24j . 

Theoretically it has been found that these two formalisms are closely related to each other, and it is shown that they 
are equivalent in the overlapping transverse momentum region where both can apply [25l427| . However, it has been 
recently observed that the experimental proton-proton data on the SSAs of the inclusive hadron production appears 
incompatible with the Sivers data from SIDIS process [28l - 30] , if one assumes that the SSAs of the inclusive hadron 
production come entirely from the Sivers contribution. This is known as the "sign mismatch" . Whether this finding 
reflects the inconsistency of our theoretical formalism is a very important question and needs to be further explored 
both theoretically and experimentally. Since the inclusive hadron production has the complication from the Collins 
contribution, the measurement for the SSAs of single inclusive jet and direct photon production in proton-proton 
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collisions [28|, |3l| could be very helpful in studying the sign mismatch, as they are free of complication from the 
fragmentation process (or the Collins effect). 

Even though direct photon production is ideal in the theoretical sense for further exploring this "sign mismatch" , 
there are no true direct photons. Direct photons and fragmentation photons are two indistinguishable contributions 
in the usual collinear factorization formalism 32], which are designated "prompt" photons. In experiments one might 
apply the photon isolation cut to reduce the fragmentation contribution, however the asymmetry measurement might 
suffer from the low photon event rates after such a cut. In any case, it is important to assess how the fragmentation 
contribution might affect the asymmetry of the prompt photons. This is the main purpose of our paper. While the 
direct photons receive only the Sivers type of contribution for the asymmetry, the fragmentation photons could receive 
both the Sivers and Collins contributions. We perform a model calculation for the quark-to-photon Collins function, 
which is then used to estimate the Collins asymmetry for fragmentation photons. 

The rest of our paper is organized as follows. In Sec. II, we give the overview on the various sources for the SSAs of 
prompt photon production. In Sec. Ill, we present our detailed model calculation for the quark-to-photon unpolarized 
fragmentation function and Collins function, and estimate their relative size. In Sec. IV, we make phenomenological 
study for the SSAs of prompt photon production by including all the sources studied in our paper. We conclude our 
paper in Sec. V. 



II. SINGLE TRANSVERSE SPIN ASYMMETRY OF PROMPT PHOTON PRODUCTION 

A. Unpolarized prompt photon production 

We consider the prompt photon production in hadronic collisions, A(Pa,s±) + B(Pb) — > j(P-y) + X. Here A is 
a transversely polarized proton with spin vector sj_, and B is an unpolarized proton. The spin-avcragcd differential 
cross section of prompt photon production contains both direct and fragmentation contributions, 

da da A " dcr frag 

e ^ = e ^ +e ^ifp;- (1) 

At leading order, the direct contribution is given by 

„ da dlr a em a s \ -» f dx' , , , v f dx „ , . „ ir ~ ~ ~ *\ 

^7^r = ^E / -jh/BW) J -f a/A (x)H^ 7 (s,t,u)S(s + t + u) , (2) 

^ a,b 

where s = (Pa + Pb) 2 , fa/A(x) an d fb/B(x') are the spin-averaged parton distribution functions, s, t, and u are the 
usual Mandelstam variables at the parton level. H^ b ^. are the well-known partonic hard-scattering functions for 
direct photon production (33l. l34j. At the leading order, they are calculated from the partonic channels qg jq and 
qq — > 7(7, and the typical Feynman diagrams are shown in Fig. [1] 





FIG. 1. Typical Feynman diagrams for direct photon production at leading order: left for qg — > 75 and right for qq — > jg. 

For fragmentation photons, in the usual collinear factorization formalism at leading order, we have 2 — >■ 2 scattering 
process to produce a parton which then fragments into a photon, with the typical Feynman diagrams shown in Fig. [3] 
The differential cross section is given by 

„ dcr frag ai \ ^ f dz _ . . f dx' „ , ,. f dx „ . . r , f - - ~ •> 

7 rf^ = ^^ / ^ Dc ^ Z) I — fb / B ^ I —fa/A(x)H^ c (s,t,u)S(§ + t + u), (3) 

where Z? c _,. 7 (z) is the quark-to-photon fragmentation function, and HK are the well-known partonic cross section 
to produce a parton [33|, [35[ . 

To see the relative contributions of direct and fragmentation photons, we define the following direct ratio 

p dljdiI 

^7 rfO^ "f" ^7 



3 




FIG. 2. Typical Feynman diagrams for fragmentation photon production at leading order. 

In Fig. 02 we plot the direct ratio R as a function of Feynman xf sX forward rapidity y — 3.5 at RHIC energy 
y/s = 200 GeV. We give the result for both the leading order and next-to- leading order calculations [32]. We find 
that the fragmentation photons actually contributes to around 50% to the total prompt photon production. Thus it 
is important to assess the effect of fragmentation photons on the asymmetry of the prompt photon production. 




FIG. 3. The direct ratio defined in Eq. ((4]) is plotted as a function of Feynman xf at y = 3.5 and y/s = 200 GeV. The solid 
line is for leading order calculation, while the dashed line is for next-to-leading order calculation. 



B. Spin-dependent cross section for prompt photon production 



In order to compute the asymmetry of prompt photon production, we need the spin-dependent cross section 
Aa(s±) = [c(s_l) — cr(— sj_)] /2, which will also contain both direct and fragmentation contributions, 



dAa „ dAcr dir „ dAcr fra s 



E. 



7 



7 7 u, ± 1 

The direct contribution contains only the Sivers type of effect, as given by [20I [26 



E. 



dA- d 
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° -e O^PP acmCts 



d 3 p 



a.b 



Ta tF {x, x) ~ X — Ta tF {x, x) 



(5) 



ab— >7 



(s, £, u)S [S + t + u) 



(6) 



where the hard-part functions H°^.^ contain the relevant initial-state interactions between the active parton and the 
remnant of the proton and have the expressions given in [20L l3(ij| . T q ^p(x, x) is the twist-three quark-gluon correlation 
function, and it is related to the quark Sivers function fyp{x, kj_) as follows [5] 



Tg iF (x,x) = - J d' 



I SIDIS, 



(7) 



where the subscript "SIDIS" here is to emphasize the Sivers function probed in SIDIS process. On the other hand, 
the spin asymmetry of fragmentation photons can receive both Sivers and Collins contributions, 



E-, 



dAcr fra S 

IppT 



1 d?P^ 



j? rfAt7 Coflins 
7 <PP^ ■ 



(8) 
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The Sivers contribution can be written as 



E-, 



dAer, 



frag 
Sivers 



£a/3 



a 



-E 



— D 



dx' 

—rh/B(x') 



dx 
x 



T a . F (x,x) 



X — T atF (x,x) 



(9) 



where H^™ 1 * represents a hard-part functions for the partonic process ab — > cd, and it incorporates both the initial 
and final-state interactions and has the expressions given in [2(| [28| . The Collins contribution for a inclusive hadron 
production has been calculated in [l3T ] , which is related to a convolution of quark transversity and quark-to-hadron 
twist-three fragmentation function. The only difference for fragmentation photons lies in the quark-to-photon twist- 
three fragmentation function, and the differential cross section is given by 



jp ^Collins 
7 tPi* 



a.b,c 



dx 



h a (x) 



dr' 



dz 

z 



,d ( H c {z) 
dz \ z 2 



1 



z x(-u) + x'(-t) 



(10) 



where h a (x) is the quark transversity, and H c {z) is the twist-three quark-to-photon fragmentation function and is 
related to the first pr-moment of the Collins function H l 9 (z,p 2 n ) 



H g (z) 



d 2 p T p T H^ q (z,p^), 



(11) 



with H^ q (z,p 2 T ) defined in the next section. The relevant hard-part function H^^ 1 ™ has been computed in [lc 
Eventually the single transverse spin asymmetry An is computed from the following definition 



i-N 



jp dAa 
rp da 



(12) 



where the spin-dependent and spin-averaged cross sections are given in Eqs. JT]) and ([5]), respectively. To calculate An 
numerically, we need the information for the twist-three quark-gluon correlation function Tq t p(x, x) and twist-three 
quark-to-photon fragmentation function H q {z). The information of T q> F(x,x) has been directly extracted from the 
proton-proton data [2(1 HU, or indirectly from the SIDIS data by using the Eq. (|7j) [H,|33|- However, the information 
of H q (z) is completely unknown. To estimate the size of H q (z) will be the main focus of the next section. 



III. QUARK TO PHOTON COLLINS FUNCTION 



In this section, we perform model calculations for photon fragmentation functions, including both the unpolarized 
fragmentation function and the Collins function. We first study the transverse momentum dependent quark-to-photon 
fragmentation functions, and then integrate over the transverse momentum to obtain the relevant unpolarized collinear 
fragmentation function D q ^.~ l {z) and the collinear twist-three fragmentation function H q (z). 



A. Transverse momentum dependent fragmentation functions 

Photon fragmentation function can be calculated from the correlation function A(z, kr) [1, H, HH, 

A(z,k T ) = ^|^e*f<0|^(Ol7X><7X|^(0)|0>| f - =0 , (13) 

where the usual gauge link is suppressed, and we have assumed that the photon is moving in — z direction with 
momentum p M = p~n^ and light-cone vector n M = [0 + ,l~,0jJ. The fragmenting quark has momentum k, with 
k~ — p~ I z and fc<r the transverse component with respect to the photon momentum p. We define pt as the photon 
transverse momentum with respect to the quark, which is related to kx as: px = —zkx- Here for our purpose we only 
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keep the terms relevant to the quark to unpolarized photon fragmentation. Then the correlation function A(z,kr) is 
given by USS, 



A(«,fcr) = - 



D, 



9^7 



{z,p 2 T )4 + Hi q (z,p 2 T )u^k T ^ 



(14) 



D q ^{z 1 p 2 ^) is the usual unpolarized quark-to-photon fragmentation function, and H 1 q (z,p 2 n ) is the quark-to-photon 
Collins function in agreement with the "Trento conventions" 13911 . For the most general case where photon's polariza- 
tion is also specified, there are more terms in the expansion [38| . We can easily project out these functions 



D q ^{z,p 2 T ) = -Tr[A(2,fc T )$], 
e^k Tu H^(z,p 2 T ) = ifr [A(z,k T )*a^n^ 5 ] , 



(15) 
(16) 



where n M = [1 + , , OjJ is a light-cone vector conjugate to n M . 




FIG. 4. The Feynman diagram which contribute to the unpolarized quark-to-photon fragmentation function D q ^ 1 {z,p%) 



In our model, the tree-level diagram describing the fragmentation of a quark into a real photon is depicted in 
Fig. [4] By contrast with the pion fragmentation calculations [40l - |4^ |. the interaction between quark and the photon 
is described by a simple point interaction with coupling ie q e^ and e q the quark fractional charge. In the actual 
calculations, we will choose light-cone gauge n ■ A cn i — for the photon field 43] to avoid photon eikonal phase 
[44j . On the other hand, we still use the covariant gauge for the gluon field, thus eikonal phase for gluon field still 
exists in our calculations. In such a set-up, we have only one Feynman diagram (at leading order) for unpolarized 
quark-to-photon fragmentation function D g _> 7 (z,p|,), as shown in Fig. 0J Thus, the photon polarization sum is given 
by 



A 



e^(p,A)e*>,A) = -<r + 



p>*h v + p"n 
n ■ p 



The calculation is straightforward, and we obtain 



D 



9^7 



(Z;Pt) 



'2ir 2 z 2 (l-z) 



l + (l-z) 2 



2zm 2 



(k 2 



v 



(17) 



(18) 



where a cm is the electro-magnetic coupling constant, m q is the quark mass, k 2 is the virtuality of the fragmenting 
quark, and it is related to photon transverse momentum px as follows: 



k 2 = 



P 2 i 



z{l-z) 



1 



(19) 



The Collins function receives contributions only from the interference between two amplitudes with different imag- 
inary parts. Since the tree-level amplitude is real, the necessary ima gina ry parts will be generated by the inclusion 
of one-loop corrections. Here we study the case of gluon loops [40h42L |45| . The relevant Feynman diagrams are given 
by Fig. El The double line in Fig. [5^c) and (d) represents the eikonalized propagator, which give rise to the factor 
l/(— n ■ t ± ie) [13, |4l|. The calculations are much more involved than the unpolarized fragmentation function, but 
nevertheless similar to those for the quark-to-pion Collins functions calculated in [4(J l4ll. In particular we note that 
the contribution from Fig. [SJd) are due to poles on the gluon and incoming quark |40l . Ell . [45[ signaling that the 
photon Collins function is universal [451. |46|. Here we give only the final results, 



Ht g (z,p 2 T ) 



>ot c 



2n 2 k 2 



H 



x(fig 



-(fig-b) 



H 



-i(fig.c) 



H 



-L(fig.d) 



(20) 




where the four terms in the bracket correspond to the four diagrams in Fig. [5] and they are given by 

,2 \ 



H 



J-(fig.a) 



2zk 2 \ k 2 r 



(21) 



H 



-L(fig-b) _ 



if 



Hi 



-L(fig.d) 



(1 — z)(/c 2 — to 2 ) 



-L(ag.c) _ 



0. 



1 



(1 - z)k 2 



m 2 ( k 2 \ 1 / m 2 mt 



(l-z)fc 2 , /(l-z)A: 2 
In ' 



(1 — z)k 2 — m 2 



(22) 

(23) 
(24) 



B. Collinear fragmentation functions 

The collinear (integrated) unpolarized fragmentation function _D 9 _s. 7 (z) is defined as 

D q ^(z)=n dp%D q ^{z,p%). (25) 

Jo 

Following [10, IH, E3' we take the upper limit p\ max to be set by a cut-off on the fragmenting quark virtuality /i 2 , 
where k 2 < /i 2 . From Eq. (|T5)) . this corresponds to 

PTmax = - Z )H 2 - zm r (26) 

Then the analytic result for D q ^ tl {z, ji 2 ) is 



D q ^(z,fi ) = e,— 



l + (l-z) 2 (l-z)(u 2 -m 2 ) / to 2 l-z 

— In ^-t; - + 2 ' 



Similarly from Eq. (|11[) . we define the twist-three fragmentation function H q {z,^ 2 ) as 



zto^ \ ju 2 — TO 2 

2\ 



(27) 



7T 







H q {z,^) = — dp 2 T p z T Ht q {z,p 2 T ). (28) 



Now let us estimate the relative size of twist-three fragmentation function H q (z,p, 2 ) compared to the unpo- 
larized fragmentation function D q ^ > . 1 {z, /i 2 ). In Fig. [51 we present numerical estimates for the analyzing power 
H q (z, /i 2 )/ D q _). 7 (z, n 2 ), separately for each of the diagrams of Fig. [5] We choose quark mass m q — 300 MeV, and 
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FIG. 6. The ratio of H q (z, fi 2 )/D 9 _> 7 (z, fi 2 ) at scale fi = 1.5 GeV is plotted as a function of z. The magenta dotted curve is 
the contribution from Fig.(5ja), the green dot-dashed for Fig.[5jb), the blue dashed for Fig.[5jd), and the red solid curve is the 
sum. 



-0.02 
-0.04 



plot this ratio at scale /z = 1.5 GeV as a function of z. The magenta dotted curve is the contribution from Fig.(5Ja), 
the green dot-dashed for Fig. [5jb), the blue dashed for Fig. [jjd), and the red solid curve is the sum. We find that 
there is a strong cancellation between the contribution of diagrams (a) and (b), similar to the quark-to-pion Collins 
function [40j ] - Thus the sum is dominantly given by the contribution from diagram (d), the gauge box diagram. We 
also notice that the quantity H q (z, [i 2 )/ D q -> 7 (z, /i 2 ) for photon case is much smaller than the same quantity for pion 
case as estimated in Ref. [40j . This leads to a much smaller Collins asymmetry for fragmentation photon production, 
as shown in the next section. 



IV. PHENOMENOLOGY 



In this section, we will estimate the SSAs of the prompt photon production in the forward rapidity region at RHIC 
energy. In order to assess the contributions from the fragmentation photons, beside the overall spin asymmetry 
An defined in Eq. (1121) . we will define the following additional asymmetries: the spin asymmetry for direct photon 
production A N r , the spin asymmetry for fragmentation photons A N &S . That is, 



E dAad " E 



A$ = -LJ%-, A%*= 7 f?'» . (29) 

a l d?P, t ^1 d?P-, 

For the fragmentation photons, there are both Sivers and Collins contributions for the spin asymmetry, we thus further 
define the Sivers asymmetry for fragmentation photons A ^ a | ivel . s , and the Collins asymmetry for fragmentation photons 



iV, Collins) 



/(frag _ 7 d3p -i (on) 

^-N, Collins _ „ J„frag ■ {° v ) 



_ 1 d^P 1 .f ra g ^7 

JV, Sivers j-, da lr ^ ' N, Collins j-, da lTl 

fra; 
JV 



Note that the spin asymmetry for fragmentation photons is the sum of Sivers and Collins asymmetry as A 
^^Sivcrs + -^jVCoiUns' However, the overall spin asymmetry for prompt photon An ^ A^ 1 ' + A N ag . 

Quark-to-photon fragmentation function has been extracted from the phenomenological study, see e.g., Ref. 48]. 
This parametrization has been used to describe the unpolarized prompt photon production at RHIC energy [49(. 
Thus, to compute the fragmentation photon cross section in spin-averaged proton-proton collisions, we will use this 
phenomenological parametrization instead of the model result in Eq. (|27p. On the other hand, since there is no 
experimental information at all for the quark-to-photon twist-three fragmentation function H q (z), we have to rely on 
our model calculation in order to estimate the Collins contribution to the asymmetry of fragmentation photons. In 
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this case, we will assume that our model calculations give a reasonable estimate on the relative size for H q (z, /i 2 ) and 



Thus we will use the following approximation, 
H q (z 



D, 



9^7 



(z,^ 2 



phenomenology 



D, 



1^1 



(31) 



odd 



and ^7v a |i vcrs , we need the twist-three quark-gluon correlation functions 



where H q (z,/j, 2 ) and D q ^(z, fj, 2 ) on the right-hand side are given by the expressions in Eqs. ([27)) and (|28l) in our 
model calculations, D q ^ 1 (z, [i 2 ) on the left-hand side is the phenomenological parametrization from Ref. [481 ] , and 
H q (z, /i 2 ) in the numerator on the left-hand side will be the quark-to-photon twist-three fragmentation function to be 
used in our calculation for the asymmetry ^^f,^^ °f fragmentation photons. For quark transversity distribution 
h a (x), we take the parametrization from Ref. [50| . 

On the other hand, to calculate 
T a j p(x, x). This function has been extracted directly from the inclusive hadron production in proton-proton collisions 
[20(, which will be labeled as "KQVY" parametrization in our plots. T a p(x, x) can also be computed indirectly from 
Eq. ([7]) with the quark Sivers function extracted from SIDIS process (36l . |37| . Such indirectly obtained parametrization 
for T q ^(x : x) from [36| will be called "old" parametrization, while that from [33] will be labeled as "new" parametriza- 
tion in our plots. It has been found in [28| that the directly and indirectly obtained T q ^p{x,x) have conflicting signs, 
for both u and d quark flavors. The future prompt photon production hopefully could help us pin down the sign and 
magnitude of T qt p(x,x). 
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FIG. 7. Single transverse spin asymmetry for prompt photon production, p' + p — > 7 + X, is plotted as a function of Feynman 
xf at rapidity y = 3.5 and center-of-mass energy yfs — 200 GeV. Left panel: the asymmetry for the fragmentation photons. 
The black solid curve is the Collins asymmetry A^^y^. Dashed curves are the Sivers asymmetry A 



frag 
1 iV,Si- 



with the red 



curve for "KQVY" parametrization, the blue curve for "new" parametrization, and the green curve for "old" parametrization 

,Co 



for T q: F(x,x). For each set, the solid curve is the asymmetry for fragmentation photons -A^ 3,8 , which is the sum of ^^coiiin 



and A^ a f ivers . Right panel: the asymmetry for the prompt photons. For each set, the dashed curve is the direct asymmetry 
Af} 1 , the dotted curve is the fragmentation asymmetry A f ™ s , and the solid curve is the overall spin asymmetry Am- 



In Fig. deleft), we plot the spin asymmetry for fragmentation photons as a function of Feynman xf at forward 
rapidity y — 3.5 and RHIC energy y/s = 200 GeV. The black solid curve is the Collins asymmetry ^47v a c ollins ■ We find 
that the Collins asymmetry is very small in the whole xf region, less than 1%. Dashed curves are the Sivers asymmetry 
^Tvfivcrs' w hh the red curve for "KQVY" parametrization, the blue curve for "new" parametrization, and the green 
curve for "old" parametrization for T q ^p(x, x). For each set, the solid curve is the asymmetry for fragmentation photons 
A^ as , which is the sum of ^jvCoiiins an< ^ -^i^Iivers- I n Fig-Height); we plot the spin asymmetry for the prompt photons. 
For each set, the dashed curve is the direct asymmetry ^4^ r , the dotted curve is the fragmentation asymmetry ^4^ as , and 
the solid curve is the overall spin asymmetry Ajq. We find that the spin asymmetry for fragmentation photons A { ™ s has 
the same sign as the direct asymmetry A^ r , thus the overall spin asymmetry An has the same sign as A^ r and ^4^ ag . 
Since the Collins asymmetry for fragmentation photons is very small, the asymmetry of prompt photon production 
can be a very good probe for the twist-three quark-gluon correlation functions T 9i p(x,x). We urge the experiments 
to measure the asymmetry of prompt photon production at RHIC. It will provide important information on the 



9 



twist-three quark-gluon correlation functions, a quantity much needed to verify our current theoretical formalism for 
describing single transverse spin asymmetry in proton-proton scatterings. The measurement can go a long way to 
resolving the so called "sign mismatch" [28M30j | . 



V. SUMMARY 



We have studied the single transverse spin asymmetry of prompt photon production in high energy proton-proton 
scattering including the contributions from both the direct and fragmentation photons. While the asymmetry for direct 
photon production receives only the Sivers type of contribution, the asymmetry for fragmentation photons receives 
both the Sivers and Collins types of contributions. In order to estimate the Collins asymmetry for fragmentation 
photons, we perform a model calculation for quark-to-photon Collins function. Within our model calculation, we find 
that the Collins asymmetry for fragmentation photons is very small in the whole kinematic region, thus the single 
transverse spin asymmetry of prompt photon production is mainly coming from the Sivers asymmetry in direct and 
fragmentation photons. We make predictions for the prompt photon spin asymmetry at RHIC energy, and find that 
the asymmetry is sizable. The asymmetry of prompt photon production can then provide a good measurement for 
the important twist-three quark-gluon correlation function, which is urgently needed in order to resolve the "sign 
mismatch" puzzle. We urge the experiments to measure the asymmetry of prompt photon production at RHIC in the 
near future. 
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